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Using lattice techniques we investigate the generation of long range cosmological magnetic fields 
during a cold electroweak transition. We will show how magnetic fields arise, during bubble 
collisions, in the form of magnetic strings. We conjecture that these magnetic strings originate 
from the alignment of magnetic dipoles associated with EW sphaleron-like configurations. We 
also discuss the early thermalisation of photons and the turbulent behaviour of the scalar fields 
after tachyonic preheating. 
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1. Introduction. 

There have been many theoretical attempts to explain the origin of large scale cosmological 
magnetic fields (LSMF) [0]. The main difficulty resides in understanding their correlation scale 
which ranges from the size of galaxies to clusters and super-clusters with an amplitude of the order 
of micro-gauss, pointing to a primordial origin. Following the work initiated in [Q], we address this 
issue in the context of a cold electroweak transition taking place after a period of hybrid inflation. 
The EW transition has been in the heart of many proposals to address magnetogenesis, linking it 
in many cases with the generation of the baryon asymmetry. The results presented here resemble 
the mechanism proposed by Vachaspati connecting the appearance of magnetic fields to that of 
sphalerons and Z-strings [Qj. 

The model we have considered is a hybrid inflation model with the bosonic field content of the 
Standard Model coupled, via the Higgs field, to a singlet inflaton: 



(D^D^] + l -(d^) 2 -V{<S>,x) (1.1) 



V(*,Z) = V + ^(g 2 X 2 -m 2 ) l^l 2 + jl^l 4 + \v 2 X 2 (1.2) 

The couplings are fixed to the standard model values for several ratios of the Higgs to W masses. 
For concreteness, we have fixed the inflaton to Higgs coupling by the relation: g 2 = IX. To solve 
the time evolution of the system, starting at the end of inflation, we have performed a numerical 
evolution based on a suitable classical approximation (more details can be found in Refs. @,|,il). 
In this work we discuss the results for Mh = 4.65 Mw- Results for more realistic values will be 
presented in [B]. 



2. Simulation results. 

The evolution of the scalars shows two regimes that determine the behaviour of the system. 
The first one is the spontaneous symmetry breaking (SSB) stage which proceeds via the generation, 
expansion and collision of bubble-like structures in the Higgs field It is in this regime where a 
rich phenomenology can arise, including the generation of the baryon asymmetry of the Universe 
and a stochastic background of gravitational waves, see Refs. [f|, ^, 0, and where we place the 
origin of the LSMF. In our simulations this regime starts right after the end of inflation (mt ~ 5 — 10 
for Mh = 4.65 Mw) and lasts until times of order mt ~ 35—40. After this relatively short stage, a 
period of slow approach to equilibrium follows where two important phenomena appear: turbulence 
and thermalisation of high momentum photons. 

2.1 First regime: Generation of magnetic fields. 

The simulations show that magnetic fields are predominantly generated in the SSB regime. We 
have already mentioned that bubble collisions are responsible for SSB. Bubbles originate at peaks 
of the initial Higgs Gaussian random field and expand from there on [Q]. As exhibited in Fig. |T] 
(Left), the regions between bubbles remain for a longer time in the false vacuum, giving rise to 
string-like structures of local minima of the Higgs-field. Strongly correlated with them, magnetic 
fields appear dominantly in the form of magnetic strings, see Fig. [I] (Right). We conjecture that 
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Figure 1: Top: (Left) Locus of points where the value of the Higgs norm is below 0.53. (Right) Locus of 
points where the U(l)-magnetic energy is above 1.5. Bottom: (Left) Smeared density of the U(l)-magnetic 
energy. (Right) Smeared density of the magnetic component of the Z field. Both above 0. 1 . All at mt = 15. 



these magnetic strings originate from sphaleron-like configurations which are copiously produced 
during SSB at the location of Higgs-field zeroes [ffl, Sh. For non-zero Weinberg angle the magnetic 



field generated by the sphaleron looks like a magnetic dipole Qili [12 ]. Then, as in a ferromagnetic 
material, these dipoles would tend to align along the string of Higgs-zeroes generating a magnetic 
string. One remark has to be done concerning Fig. |]. The plot has a lower cutoff giving the 
impression that magnetic lines are open. This is just a projection effect, and places where the 
magnetic strings seem to end are just places where the magnetic flux spreads. 

The picture that we extract resembles the mechanism proposed in Ref. [§], in which magnetic 
strings also appear at the electroweak phase transition, and are related to Z-strings. The bottom 
half of Fig. [I] shows the smeared Z-magnetic field together with the smeared U(l)-magnetic field 
at mt = 15. The correlation is manifest. In this plot the string structure has been enhanced by 
performing smearing (averaging the fields over the neighborhood of a point). This helps in two 
ways: first, it smooths the central singularities of the strings and secondly, it clears up the radiation 
background of non-coherent fields. The procedure respects the regions where the scales are roughly 
constant and does not disturb the whole picture as far as the neighborhood involved in the averaging 
is smaller than the width of the string (~ 1 /Mw). 

Having a mechanism for generating magnetic fields, two important questions arise, and the 
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Figure 2: Left: Averaged Bi(ro), Eq. (2. 1), at mt = 10. The center of the box is chosen to be the maximum 
value in the center of a string. The plateau length provides an estimate of the minimal length of a representa- 
tive string. Right: Dependence of B\ = (Bj) on the time-interval length T, see Eq. (2.2). B\ is normalized 
to the total energy density. 

viability of this mechanism as LSMF generator depends on their answer. These questions are: 

• Are the scales of the generated magnetic fields large enough to be the seed of LSMF? 

• Do they persist in time long enough to be used as seeds of LSMF? 

A qualitative answer to the first question can be extracted from Fig. |]. Individual strings 
merge due to magnetic reconnection forming a cluster of the order of the box size. This merging 
of individual strings is essential to enhance the large scale structure of the fields. In order to give a 
quantitative measure of the nature of the magnetic field objects as actual strings, we have analysed 
the following quantity: 

B L (r Q ) = ^[ dx 3 B 2 (x), (2.1) 
£ JL(r ) 

where the integration is on a box of length L, centered at a point r belonging to one of the strings. 
Figure |2| (Left) shows the L dependence of L 2 Bi(ro), averaged over several configurations. The 
expected behaviour of Bi{ro) if ro is in the central core of a string-like object is clear. Inside the 
string the magnetic intensity is roughly constant, so the integral goes like the volume of the box and 
Bt( r o) is constant. Once the size of the box gets larger than the diameter of the string, the integral 
only grows along the direction of the string and Bi(ro) ~ L~ 2 . This gives the plateau seen in the 
figure. The small errors indicate that the string objects are quite configuration independent.The 
length of the plateau provides an estimate of the minimal length of the string as an isolated object. 
It indicates a mean minimal string length of order M^L ~ 4 — 5, which is a significant part of the 
total length. For larger values of L the merging of strings and the fact that more strings are entering 
the box, change the behaviour to a volume like one, with 5l(to) roughly constant. 

The time persistence analysis is quite more subtle. A detailed analysis of the behaviour in 
time of the low momentum magnetic field spectra is on the way and will be presented in [QJ. Some 
preliminary conclusions can already be presented here. A time average can be defined as follows: 

B T (to) = 4 / dtB(t), (2.2) 

T jT(t ) 
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Figure 3: Scalar variances, Eq. (2.4), plotted in logarithmic scale for different lattice spacings and volumes. 



this is, the integral over a time interval of length T that starts at to- This time-averaging can be 
seen as a way to eliminate short scale magnetic fields, as the radiation fields. For instance, for a 
monochromatic wave, 

(\B T ?) ~ I J dk\B k ? W\kJ) , W = SiH ^ 2 \ (2-3) 

where the average is both over space and configurations. For large T, momenta with \k\ » (T)~ l 
are suppressed. We can use that to estimate the length of the time interval for which all momenta 
but our minimal momenta p m m = 2n/L are suppressed. This gives a window of possible intensities 
of the magnetic seeds coherent at scales larger or equal to the scale given by the size of the box. 
Fig. H| (Right) shows the dependence of (|B7-| 2 ), normalized to the total energy density (po), on the 
length of the time interval T, for p m [ n = 0.15m. Note that the average turns out to be quite indepen- 
dent of the starting point of the averaging to, indicating that the sub-leading large scales are quite 
independent of time. From the values above mT ~ 6 we obtain a window (Bj) ~ [10~ 2 , 10~ 3 ] po- 
Assuming the magnetic field expands as radiation, this would give magnetic fields today of order 
[1 , 0.3];UG, which are in the range of the observed ones in galaxies, and even in clusters of galaxies, 
where no-extra amplification through a dynamo mechanisms is expected. 

2.2 Second regime: Turbulence and photon thermalisation. 

Once the SSB takes place, the system starts a slow approach to equilibrium. This proceeds via 
a period in which the scalars and also some of the gauge degrees of freedom undergo a turbulent 
regime. Turbulence in the scalar fields manifests through the time evolution of their variances 
shown in Fig. ||[ As predicted by JTo| ] and also shown in |9|], variances follow the law: 

for m-particle interactions and / denoting each of the scalars. The change in model parameters 
with respect to [Q], where M w = 6 GeV, results in a modification of the turbulent behaviour. The 
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Figure 4: Left: log(n^ + 1) and Right: log(D(B)/B 2 ), Eq. (2.5), for the spectrum of photon radiation at 
times mt = 85 and mt = 225 respectively. 

parameter m, now turns out to be m = 5(4) for the Higgs (inflaton), which differs from the former 
m = 4(3) obtained in [Q], but is in agreement with [BJ. 

In Ref. [Q] we showed the self-similarity of the momentum spectra of the SU(2) degrees of 
freedom long after SSB. This additional sign of turbulence indicated a departure from thermal 
equilibrium. High momentum photons, however, thermalise very early in the evolution. In thermal 
equilibrium the radiation occupation number is expected to behave according to the Bose-Einstein 
distribution: rik = (exp(j3&) — 1) _1 , from where the thermalisation temperature can be extracted. 
Fig. H (Left) represents ln(l +n^ 1 ) for photon radiation. Photons with momenta above m follow 
the Bose-Einstein law. The temperature extracted is of the order of Mw from times mt ~ 40 on. 
Thermalisation is corroborated by analysing the distribution of norms of the magnetic field. In 
thermal equilibrium this distribution is given by [Q] : 

D(B) = B 2 e- 3B2 /(2<* 2 » ; (B 2 ) = ^ (2.5) 

Fig. H right shows D{B). From time mt ~ 40 on the equilibrium distribution is reached, with 
a temperature in agreement with the one extracted from the spectra. In both plots a systematic 
deviation from the thermal behaviour is observed at the low momenta, large norm, part of the fits. 
It is there where the interesting magnetic string structures reside. An analysis of the possible, late 
time, turbulent and helical behaviour of the low momentum part of the magnetic field spectrum is 
on the way and will be reported elsewhere. 

3. Conclusions. 

We have analysed numerically the proposal that long range magnetic fields could be generated 
during a cold electroweak transition after a period of low scale hybrid inflation. The generation 
mechanism is mainly based on two facts: 

• At the SSB stage bubble-like structures, associated to local maxima in the Higgs-field norm, 
appear. Points outside the bubble front, remaining close to the false vacuum, form string like 
structures. 
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• Bubble collisions give rise to sphaleron-like configurations attached to the location of zeroes 
of the Higgs field. For Ow ^ these sphalerons behave as magnetic dipoles. 

These two ingredients together lead to an allignment of the sphalerons' dipoles forming mag- 
netic string networks as we have observed in our simulations. Some results concerning the coher- 
ence and intensity of the generated magnetic fields have been presented. A further analysis of the 
time evolution of the magnetic field, as well as the study of the dependence of the Mu to M\y ratio 
will be presented in 

We have also discussed some features of the late time behaviour of the system, among them, 
thermalisation of photon radiation and turbulence in the scalar fields. 
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